Phenomenological framework for nonequilibrium steady states of molecular motors 
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An expression for the energetic efficiency of a molecular motor is presented in terms of an effective 
temperature, which was defined based on the ratio of the correlation function to the susceptibility 
of its velocity. We also present a numerical result regarding this temperature for a thermal ratchet 
model. Furthermore, our expression of efficiency using known values for kinesin can adequately ex- 
plain the experimental results. The framework presented here has a closed form phenomenologically, 
and is independent of any detailed models. 
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What is the essential difference between our muscle 
and the Carnot cycle as an energy transducer? The re- 
cent development of single-molecule-detection techniques 
Q for molecular motors led to new perspectives for this 
question, by directly visualizing their noisy way of chcmo- 
mechanical energy transduction at the single-molecule 
level. Many theoretical studies have been made based 
on such experiments to present various models Q, upon 
which energetics have also been discussed However, 
since most of them are bottom-up type theories, they 
are in a dilemma that they cannot be compared with 
experiments until a model and parameters for a molec- 
ular motor have been determined completely. This is 
the very difference from a top-down type theory, namely 
a phenomenology such as thermodynamics, which deals 
with energetics based on macroscopic observables alone. 
It thus sounds meaningful to construct a phenomenology 
for a molecular motor. 

In this Letter, we present an expression for the energy 
efficiency of a molecular motor based on the susceptibil- 
ity and the correlation function of its velocity, both of 
which are observables. For this purpose, we introduce 
an effective temperature as the ratio of the correlation 
function to the susceptibility. According to Fluctuation 
Dissipation Theorem (FDT) 0, this effective tempera- 
ture coincides the temperature of the heat bath when the 
motor is in equilibrium, while we show that it varies in 
general with respect to the timescale of measurement, if 
the motor is at a nonequilibrium steady state. 

First, we define the effective temperature. The velocity 
v(t) of a colloidal particle in solution responds linearly to 
a small external force f{t) as 
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x(t-t')f(t')dt', 



(1) 



where (...) is the statistical average, Vq is the mean veloc- 
ity of a reference system (a steady system without per- 
turbation), and x(t) is called a response function. The 
Fourier-Laplace transformation of given as 



is called the susceptibility. The 1st type of FDT states 
that if the reference system is in equilibrium, the real part 
of the susceptibility x'[H = R- e xM is proportional to the 
correlation function of the fluctuation of the reference 
system. 



$[u]=T x'[ 



(3) 



where 



1 r°° 

$H = - J (KO) - v )(v(t) - v ))e^dt. (4) 

To is the temperature of the reference equilibrium system. 
We set the Boltzmann constant to unity. 

We proceed to the case where the reference system is 
not in equilibrium but rather in a nonequilibrium steady 
state. In this case, since FDT may not hold in a normal 
sense, the coefficient in Eq. (|2J| differs from usual tem- 
perature in general even when the system is in contact 
with a single heat bath. We can adopt Eq. J3J) as a defini- 
tion of temperature for the nonequilibrium steady state. 
After calculating the susceptibility (Eq. J2J) and the cor- 
relation function (Eq. ijljl) for a nonequilibrium steady 
state as a reference system, the effective temperature is 
obtained as 



7R 



(5) 



XH = / X(t)e^dt, 



(2) 



We call this a 'dissipation temperature'. This formula- 
tion for a degree of freedom asymmetric with respect to 
a time reversal, such as velocity, is a counterpart of that 
introduced in spin glass systems for a symmetrical de- 
gree of freedom such as magnetization. It may be obvious 
that the dissipation temperature varies depending on the 
timescale uj of observation in general. 

Before an energetics argument, we compute the dis- 
sipation temperature for a typical model of thermal 
ratchet, known as a flashing ratchet 6, 7], in a nonequi- 
librium steady state, to gain insight into this newly de- 
fined temperature. A flashing ratchet is one of the well 
known models for molecular motors || , in which Brown- 
ian motion of a particle can be rectified in one direction at 
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a nonequilibrium steady state. In this paper, we suppose 
that a particle has two chemical states, each of which 
provides a periodic potential profile Ui(x) and U 2 {x), re- 
spectively (see the inset in Fig. 1 (a)). The transition 
between these two states is assumed to be stochastic at 
the same and constant transition rate a. The Langevin 
equation and the transition dynamics for the particle are 



mx + 7± = —d x Ui(x) - 

m) = o, ww)) = 

[1] 5 [2]. 
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2jT S(t-t'), 



(6) 



We assume that the mass m and viscosity 7 of the par- 
ticle are the same for each state. £(t) is white Gaussian 
thermal noise with a magnitude given by the temperature 
of the heat bath To. 

We numerically solved the equations JHJ with some ini- 
tial conditions. After some time lag, the system reaches 
a steady state. First, the double-sided power spectrum 
density is calculated by Fourier transformation for the 
velocity v(t) — x(t) of the particle under no pertur- 
bation, i.e. / = 0. By the Wiener-Khinchine theo- 
rem this amount is equal to 3>[w]. Next, the ve- 
locity is determined in the presence of a small oscillat- 
ing force f(t) = ecos(cut). The average linear power 
loss f(t)(v(t) — «o) results in the in- phase susceptibility 
x'[w] = a(f(t)(v(t) — v ))/e 2 , where a = 1 for ui = and 
a = 1/2 for cj 7^ (see Fig. 1 (a)). Finally, we can obtain 
the dissipation temperature by Eq. (jSJ). We confirmed 
that the dissipation temperature calculated in this pro- 
cedure equals To in the case of a = 0, i.e. in equilibrium. 
Figure 1 (b) shows the numerically computed dissipation 
temperature T[ui] in a nonequilibrium case. As shown, 
the dissipation temperature has different values with a 
long time limit (u> — > 0) and a short time limit (u) — > 00), 
and undergoes a sigmoid-like transition at an interme- 
diate frequency u>*. The dissipation temperature equals 
the temperature of the heat bath in the short time limit 
and increases with a long time limit |9(- The transition 
timescale 2-k/uj* coincides the mean passage time r of 
the particle in this ratchet system, i.e. the mean time to 
pass along a single period of the potential. 

We now address the energetics in terms of this dis- 
sipation temperature. We assume that the motion of 
a ratchet system or a molecular motor can be approxi- 
mated by the following generalized Langevin equation: 



r(t-t')v(t')dt' = -/ + £(t)+5(t). (7) 



— / is a small external load. £(t) is a thermal noise which 
satisfies normal FDT as = and (£(*)£(*')) = 

ToT(t — t'). Suppose that S(i) is the nonequilibrium por- 
tion of noise with the mean value (S(t)} = So and the 
correlation ((H(Q) - S )(H(i) - S )} = M(t - t'). We 
assume no correlation between and S(t). 
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FIG. 1: Numerical result for the dissipation-temperature com- 
putation for a flashing ratchet model (color), (a) In-phase sus- 
ceptibility x'M an d two-sided power spectrum density $[o;]. 
The black and red lines indicate the susceptibility and the 
power spectrum, respectively. Inset : Scheme for a flashing 
ratchet. The parameters used here were Uq = 5, xo — 0.2, 
a — 0.5, mass of the particle m = 0.01, viscous coefficient 
7=1, and the heat bath temperature To = 1. (b) The dis- 
sipation temperature T[u>] calculated from the data in (a). 
The horizontal green line represents the temperature of the 
heat bath To. The vertical red line is the inverse of the mean 
passage time r = 65.0 ± 0.874, which was also computed nu- 
merically. 



By solving Eq. J7J with Fourier transformation and 
comparing the results with the definition of x[oj\, we ob- 
tain 



xH- 1 = 



POO 

T[u>] = / T(t)e wt dt, 
Jo 



(8) 



where T[lo] is a Fourier-Laplace transformation of the 
friction function T(t). Similarly, the following relation 
is obtained, 

$M = (T r'M + AfM)/|rM| 2 , (9) 

where M[w] = 1/2 f™^ M(t)e iuJt dt, and T'[uj} = ReF[w]. 
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From Eq. © (note that x'M = r>]/|I>] | 2 ), 

MM 



T[u] = T 



(10) 



This relation corresponds to the 2nd kind of FDT. We 
define ATM =T[u]-T = M[uS\/T'[uf\. 

We now perform long-time averaging on Eq. 10 after 
multiplying by v(t) and obtain the expression 



/M*)> = - 



i>] dto 

-co r[w] °2tt 



r(t-i')M*M*'))d*' 

E (v(t)) 



(11) 

do; 



This expression represents a balance of energy flux. The 
left-hand side represents the work against the external 
load. On the right-hand side, the 1st term corresponds 
to dissipation into the heat bath, the 2nd term is the 
energy input from the bath, and the 3rd and 4th terms 
reflect external energy input. Therefore, the efficiency of 
energy transduction from external energy to work can be 
defined as 



f(v) 



(12) 



where Q ir = {T'[u]/T[tJ\) AT[w]dw/27r. It is clear 
from this expression that a dissipation temperature 
greater than the heat bath temperature leads to irre- 
versible dissipation and decreases the efficiency. On the 
other hand, in the ideal case where ATM, = and 
/ — > Hq, which corresponds to a quasistatic process, we 
get e — ► 1. Thus, Eq. i(T2"|> represents how much work 
can be extracted from thermodynamical free energy |l0| . 
This definition of efficiency must be identified from what 
usually appears in the theory for the Carnot cycle in ther- 
modynamics, rj, which reflects how much work can be ex- 
tracted from an absolute amount of input energy. In the 
present context, 77 is expressed as 77 = f(v)/ (Q + E (v}), 
where Q = (r'M/rMJ T[u]du/2ir. 

It is also clear that the condition Q u > is necessary 
because the efficiency e must be less than unity. Thus, 
the fact that the dissipation temperature is equal to or 
greater than the heat bath temperature is associated with 
the stability of the nonequilibrium steady state. 

Let us compare this result with the experimental re- 
sults for the conventional kinesin. Although neither the 
correlation function nor susceptibility have been mea- 
sured for kinesin, information is available regarding the 
diffusion coefficient, mobility and so on. Thus, we adopt 
the following approximation based on the assumption 
that the dissipation temperature of kinesin has the same 
characteristics as that of a flashing ratchet (Fig. 2 (b)). 
First, we suppose that the dissipation temperature of ki- 
nesin is also a sigmoidal function with a transition around 
the mean passage time r. We also suppose that the 



TABLE I: Physical parameters of kinesin under a saturated 
ATP concentration (> 1 mM) at room temperature [Til Il2j|. 
fi, D, t, and Ho represent the mobility, the diffusion coef- 
ficient, the mean passage time, and the internal force, re- 
spectively. In these experiments, a small polymer particle 
attached to kinesin molecules was trapped with optical tweez- 
ers and the movement of the particle on a microtubule was 
observed at nano-meter and milli-second resolution pH . 



fi (nm/s-pN) 


D (nm 2 /s) 


r (ms) 


H ( P N) 


100 


1300 


10 


8 




FIG. 2: Energetic efficiency against a load force for kinesin. 
The solid line shows the approximation given by equation (|13|l 
with known physical values for kinesin. Black points are the 
experimental data from Ref. [T^]. 



imaginary part of T[u>], which is related to the inertia 
of the motor, is negligible above the timescale of inertia 
m/7, which is much smaller than r. Under these assump- 
tions, Q lr can be approximated as Q„ s» 2AT[0]uj* /2ir = 
2(D/n — To)/t, where D and (x are the diffusion coeffi- 
cient and mobility, respectively. Thus, efficiency can be 
expressed as 



2(T>//x-To)/r + M Ho(Ho-/)' 



(13) 



We used the values shown in Table 1 from Refs. 
[ill . The value of H was determined from a stalling 
force for kinesin, since the force- velocity curve of kinesin 
is almost linear. Based on these values, the long-time dis- 
sipation temperature of kinesin is T[0] = D//cbM ~ 900 
K, while T) ~ 300 K (We only get back Boltzmann con- 
stant /cb in this paragraph). Equation (|13f) can be evalu- 
ated as e w /(8 - /)/ (16.8 + 8 • (8 - /)). Figure 2 shows 
this function along with the experimental data for kinesin 
taken from Ref. [12j . These data agree very well without 
any fitting parameters. 
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A remarkable feature of the present theory is its closed 
formalism only with observables. We assumed no specific 
model to calculate the efficiency. However we need to 
discuss the assumption used here: the generalized phc- 
nomenological Langevin equation (JJJ. First, we must 
carefully consider whether the motion of a molecular mo- 
tor is a Markovian process that can be described with 
such a linear Langevin equation. This assumption re- 
quires that the internal force of the motor H(t) must be 
small enough. For the same reason, the present theory is 
limited on the case of a small external load — / 01 . This 
assumption seems to be satisfied for kinesin, since its re- 
sponse to a external force is linear. Moreover, although 
we assumed an equation only for velocity, it is neces- 
sary that there is little or no dissipation for the other 
degrees of freedom, so that the above assumption docs 
not affect the energetics. If this point is violated, an 
argument for another degree of freedom is conducted in 
parallel with the argument presented here, by measuring 
the dissipation temperature for such a degree of freedom 
and constructing a suitable generalized Langevin equa- 
tion. In this case, the dissipation temperature will be a 
tensor (we do not present the details here). However, the 
fact that the present theory which considers only veloc- 
ity can reproduce the experimental data of kinesin might 
mean that dissipation for the other degrees of freedom is 
negligible compared to that for velocity. 

In conclusion, we have presented an expression for the 
energetic efficiency of a molecular motor using observable 
quantities alone. A concept of the dissipation tempera- 
ture is also introduced based on the susceptibility and the 
correlation function of the motor. The dissipation tem- 
perature for a flashing ratchet exhibits two regimes with 
respect to the observation timescale. In the short-time 
regime, the dissipation temperature equals the tempera- 
ture of the heat bath, while it increases in the long-time 
regime. The two regimes meet around the mean passage 
time of the ratchet system. The expression of efficiency 
derived in the present paper showed a good agreement 
with the experimental data for kinesin. 

The dissipation temperature is an observable quan- 
tity in principle using recently developed single-molcculc- 
detection techniques for molecular motors, although 
there has been no previous observation. The measure- 
ment of the dissipation temperature of processive motors, 
such as kinesin, or cooperative systems, such as muscle 
fibers or flagella, should be an interesting problem in it- 
self, to clarify their energetic properties over a wide range 
of timescale. The data will be also useful for testing pre- 
viously described models of molecular motors. 

The issue of effective temperature is also related to fun- 
damental problems of thermodynamics in noncquilibrium 
systems. It has been recently reported that the effective 
temperature based on FDT actually has the nature of 



temperature in several nonequilibrium systems |14l Il5| . 
We need a further investigation to judge whether the dis- 
sipation temperature discussed here can also be a natural 
temperature for molecular motors, which hopefully leads 
us to an understandng of nonequilibrium thermodynam- 
ics for molecular motors. 
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